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THE  TRANSVERSE  RESPONSE  OF  THE  LUMBAR 
SPINE  UNDER  LONGITUDINAL  LOADS 

H.  E.  Krause  and  M.  Shirazi 

University  of  Dayton  Research  Institute 
Dayton,  Ohio 


ABSTRACT 


A  novel  continuous  model  of  the  spine  is  presented.  The  trans¬ 
verse  motion  in  the  sagittal  plane  of  the  spine  of  sitting  human  subjects 
exposed  to  vertical  vibrations  revealed  considerable  bending  along  the 
lumbar  spine  and  negligible  bending  along  the  thoracic  spine.  There¬ 
fore,  the  model  consists  of  a  curved  rod,  representing  the  lumbar  spine, 
longitudinally  loaded  by  a  mass,  representing  the  thorax.  The  differen¬ 
tial  equation  of  the  transverse  motion  was  derived  and  solved  by  making 
a  product  assumption.  No  transverse  displacements  and  no  bending 
moment  was  assumed  at  the  pelvic  end.  At  the  thoracic  end,  a  shear 
force  and  a  bending  moment  are  applied,  representing  translatory  and 
rotatory  inertia  of  the  attached  rib  cage. 

Eigenfunctions  and  eigenvalues  depend  on  longitudinal  loading.  The 
eigenfunctions  correlate  well  with  data  obtained  over  a  large  range  of 
experimental  conditions.  The  solution  for  the  time  distribution  contains 
various  distinct  harmonic  components  if  an  external  force,  alternating 
sinusoidally  at  only  one  discrete  frequency,  is  applied.  This  effect  is 
due  to  a  periodic  parameter  in  the  differential  equation. 

INTRODUCTION 

Most  observed  spinal  injuries  result  from  external  forces  applied 
in  the  longitudinal  direction.  Long-time  exposure  to  moderate  external 
loads  can  develop  slowly  increasing  damages.  Compression  of  the 
spinal  column  is  the  prevailing  stress  mode. 


I  1 

I  In  recent  tests,  Vulcan  et  al.  observed  considerable  bending  of  the 

jj  spine  of  human  cadavers  that  were  exposed  to  spinal -direction  impact. 

\  Considerable  bending  of  the  lumbar  spine  of  sitting  human  subjects  under 

\  .  *  .  2 

l  vertical  vibrations  was  observed  by  Krause.  In  view  of  these  findings 

:  it  se^ms  that  transverse  displacements  as  well  as  longitudinal  displace- 

1  ments  should  be  considered. 

Ber-ding  relieves  stress  in  certain  areas  and  increases  it  in  others, 
if  it  is  superimposed  on  existing  compression.  Therefore,  bending  may 
be  the  factor  that  determines  the  location  of  injury  if  it  is  present  to  a 
significant  magnitude. 

This  is  an  exploratory  study  in  which  only  transverse  responses 
under  longitudinal  loads  will  be  considered.  In  particular,  we  attempt 
to  account  quantitatively  for  the  effect  of  a  longitudinal  load  and  initial 
curvature  on  the  bending  stress,  its  distribution  along  the  lumbar  spine, 
and  the  conditions  of  dynamic  instability.  Longitudinal  responses  will 

3 

be  neglected.  It  has  been  shown  by  Bolotin  that  transverse  vibrations 
are  significantly  disturbed  by  longitudinal  vibrations  only  at  resonance 
of  the  latter.  Longitudinal  vibrations  can  therefore  be  neglected  if  we 
assume  their  resonant  frequencies  and  those  of  the  transverse  vibrations 
do  not  coincide. 

SPINAL  MODEL  CONFIGURATION 

|  For  a  hypothesis,  a  spine  model  is  assumed  that  consists  of  two 

parts  (Figure  1).  They  are  a  curved  rod  to  represent  the  lumbar  spine, 
and  a  rigid  mass  connected  with  it  that  represents  the  thorax. 

Experimental  observations  have  led  to  this  configuration  (Figure  2). 
The  amplitude  of  the  first  derivative  of  the  transverse  displacement  of 
the  thoracic  spine  above  the  10th  thoracic  vertebra  of  a  sitting  human 
subject  under  vertical  vibrations  from  10  to  40  Hz  is  almost  constant. 

The  second  derivative  is,  therefore,  almost  zero  and  hardly  any  bend¬ 
ing  exists.  It  is  assumed  that  the  rib  cage  increases  the  bending  stiff¬ 
ness  of  that  part  of  the  spinal  column. 
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Figure  1,  Model  Configuration  for  Lumbar  Spine  and  Thorax, 


Transverse  and  Torsional  Vibration  Amplitude 
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It  should  be  noted  that  Figure  2  includes  the  torsional  deflection 
under  the  same  experimental  conditions.  A  nodal  point  can  be  observed 
around  the  6th  and  5th  thoracic  vertebra.  This  point  was  observed  for 
13  frequencies  within  the  test  range  from  10  and  40  Hz  and  for  30  ex¬ 
periments  at  each  frequency.  This  nodal  point  coincides  with  a  peak  at 
that  point  of  almost  all  published  curves  of  injury  incidence  along  the 
spinal  column. 


EQUATION  OF  MOTION 


The  equation  of  transverse  motion  of  a  curved  beam  has  been  pre- 


2u  (j2u- 

«  -  P(t)  — ys- 


where 

E  =  modulus  of  elasticity 

I  =  aereal  moment  of  inertia 

v(x,  t)  =  transverse  displacement  from  initial  curvature 

v  (x)  =  initial  curvature 
o 

P(t)  =  longitudinal  load  (compression  positive) 
p  =  mass  per  unit  length  of  rod. 

The  function  vQ(x)  =  0  if  the  beam  is  straight.  The  nonhomogeneous 
equation  (1)  reduces  to  a  homogeneous  equation 


d4v 


<3zv 


El  —  +  Pit) 
dx4  dx2 


a;2 


s  0 


A  product  assumption 

00 

vlx.t)*  Ivn(x)Tnlt) 
n*l  n 

where 

Vn(x)  =  spatial  distribution 
Tn(t)  =  time  distribution, 


(2) 


(3) 
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;  decomposes  the  partial  differential  equation  into  a  system  of  ordinary 

j  ones.  They  are 


£1 


(4) 

Vx)  + 


^V, 


n 


(4) 


Tn+aJTnn0 

P(t)  is  here  considered  to  be  constant. 
The  solutions  for  these  equations  are 
A.  Spatial  Distribution 


(5) 


Vnl*)«Din  sinr2x  +  D2n  cosr2  x 
4-  Dgn  slnhrj  x*H  04ll  coshr,  x 

The  rj  and  r2  are 


B.  Time  Distribution 


(6) 


Tn(t)  =  Aln  cosnnt  +  Atosinnnt 


(8) 
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where 

&n  =  natural  frequencies 
Dmn  =  constants 
rt ,  r2  =  eigenvalues 


The  constants  Dmn  will  have  to  be  determined  by  the  boundary  con¬ 
ditions.  Two  different  sets  of  boundary  conditions  will  be  considered  in 
the  following  sections.  A  solution  for  the  homogeneous  (that  is,  straight 
beam)  equation  of  motion  and  for  the  nonhomogeneous  (that  is,  curved 
beam)  equation  of  motion  will  be  derived  for  each  set  of  boundary 
conditions. 


RESTRAINED  THORAX 

The  thorax  is  supposed  to  be  restrained  so  that  no  transverse  dis¬ 
placement  or  rotation  can  occur.  The  idealized  model  for  the  lumbar 

spine  under  these  conditions  is  depicted  in  Figure  3.  This  configuration 

3 

has  been  treated  extensively  by  Bolotin  and  other  researchers  in  the 
field  of  dynamic  instability  with  the  inclusion  of  nonlinear  damping  and 
nonlinear  inertia. 


A.  Straight  Rod 


The  solutions  for  the  constants  in  the  spatial  distribution  equation 
for  the  straight  beam  (equation  6)  are  determined  by  the  boundary  con¬ 
ditions  which  are 


V(Z)  =  V(0)  =  0 


dx  x-t 


-  <^> 

dx  x=0 


=  0 


(9) 


Transverse  displacements  and  bending  moments  are  not  admitted 
at  either  end  of  the  lumbar  spine.  This  furnishes  the  equations: 
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V(!)  =  0  =  Dj  sinr2f  +  D2 cos  r2  i 


+  Dj  sin  h  rx  1  +  D4  cos  h  rx  1 

(10) 

V(0)  =  0 

=  D2  +  D4 

(11) 

v"(i)  =  0 

=  “  Dj  r22  sin  r2  i  -  D2  r22  cos  r2  £ 

(12) 

+  Dj  rx 2  sin  h  rt  i  +  D4  rx  2  cos  h  rx  f 


V"(0)  =  0  =  -D2r22+  D4  rx2  .  (13) 

These  are  four  homogeneous  equations  for  the  four  unknowns.  In 
order  to  obtain  solutions,  the  characteristic  determinant  has  to  be  zero. 


t  ~  0  = 


0 

0 


sin  r zl 

2  * 

rZ  smV 


cos  r .,£ 

2  . 
r2  cos  r^i 


0 

0 


sin  hr^  £ 
+r^2  sin  hr^l 


cos  hr^  £ 
Tj2  cos  hr 


2  2  2 

1  *  +  (tj  +  )  sin  r^£  sin  hrj  £  =  0 


(H) 


This  condition  is  met  if 


or 


rji  •  nn 


r 


1 


£  = 


0 


(15) 

(16) 


The  latter  condition  requires  P  =  0  and  is,  therefore,  trivial.  The 
eigenvalues  are: 
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(1?) 


.  -  £H.--v/p.  Sif.  p  ,2  .  unz 

2  l  ~Vz£i  *  V  (2EI*  +  EI 

This  equation  can  be  solved  for  the  natural  frequencies 

(2EI  (~)2  -  P)2  -  P2 
QC  -  —  j  - 

n  4EIu  *  (18) 

It  is  apparent  that  these  frequencies  depend  on  P.  If  P  is  aero,  we 
obtain 


.  it  .4 

> 


EI 

U 


(19) 


The  frequency  equation  can  be  rearranged  to  assume  the  following 
form: 


,n»M  JSL 
1  1  *  n 


&L\Z  l 

<T>  S 


(20) 


The  natural  frequencies  equal  those  of  the 

* . 


straight  beam  except  for 


The  relationship  between  Qn  and  P  is  ^rabolic.  The  e  lution  vfx.t) 
is  of  a  periodic  nature  if  Q  is  real  and  a  is  of  a  non-periodic  nature 
if  Q  is  imaginary.  The  transition  between  both  regimes  occurs  when 
Q  «  0.  This  is  the  case  if 


P  e  p 


crit 


s  a 


<|V  El . 


(21) 


For  P  >  P^rtt  we  have  buckling  and  for  P  <  P  the  beam  is  stable. 
Pcrit  ia  equal  to  the  Euler  buckling  load  for  this  type  of  beam  support. 
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The  dependence  of  the  natural  frequencies  on  the  longitudinal  load 
may  be  of  practical  consequence.  For  example,  the  response  of  the 
spine  to  external  forces  seems  to  be  different  under  different  sustained 
vertical  loads  as  they  occur  in  a  dive,  pullout,  or  tight  curve.  The 
amplitude  of  free  vibrations  is  expected  to  change  with  the  natural  fre¬ 
quency  if  energy  is  conserved  during  transition  from  one  longitudinal 
load  to  another. 

The  eigenfunctions  for  this  model  are 

Vn(x)  a  sin  r2n  x  .  (21) 

They  are  identical  with  the  eigenf tuitions  of  the  unloaded  straight  beam. 
It  is  easy  to  see  that  these  functions  are  orthogonal. 

B.  Curved  Rod 

The  initial  curvature  will  be  considered  in  the  following  paragraph 
by  solving  the  nonhomogeneous  equation  of  motion,  equation  (1). 

The  boundary  conditions  remain  the  same. 

Again,  a  product  aesumption  is  made: 

’*'•*»  *  J,  V*>  V  •  <“» 


V#|x)  represents  the  eigenfunctions  that  are  now  known  t*om  the  pre¬ 
ceding  paragraph.  Tfl(t)  are  solutions  that  we  seek.  Considering  the 
orthogonality  of  the  eigenfunctions,  this  assumption  leads  so 
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.ntt.2  P, 
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f  P  V  (X)  dx 

0  dx*  » 


I*/  vjuidx 
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(14) 
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The  initial  curvature  may  be  represented  by  a  Fourier  sine  series. 

This  series  contains  only  sine  functions  because  the  initial  curvature  can 
be  an  odd  function  with 

V  (0)  =  V  (f)  =  0. 
o  o 

The  second  derivative  of  the  initial  curvature  is 


d2v0U) 

dx2 


®  n* 

-  I  Son  s,n  7* 
n=i  * 


(25) 


which  results  in  a  set  of  differential  equations 


2  pit) i  p<n s 


t.+t.  [(•?■>  |L«a=»  ^}- 


on 


(26) 


These  equations  are  uncoupled  due  to  the  orthogonality  of  the  eigen¬ 
functions  and  due  to  the  orthogonality  between  the  eigenfunctions  and  the 
Fourier  components  of  the  initial  curvature.  With  the  initial  conditions 


dv 


vtx,o)«f.U)  ;  <■—)  M2(x)  , 

1  01  t=0 


(27) 


we  obtain  the  solution 


Tn  (»1*  jjj  (/»2(x)sin  yxd*^sinan* 

+  f—  /f.Ulsin  5p*,iX--^»a.leo5ant  + 
**  1  o  *  M  J 
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on 


(28) 
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The  longitudinal  load  P,  as  well  as  the  modulus  of  the  Fourier  compo¬ 
nents,  is  contained  in  the  coefficient  of  cos  ft  nt.  Thus  both  of  them  have 
an  effect  on  the  amplitude  of  the  free  motion. 

The  coefficient  of  the  cosine  function  of  the  n^  eigenmode  contains 

th 

only  the  Fourier  coefficient  of  the  n  component  of  the  initial  curvature. 

If  some  of  these  components  are  zero,  then  the  equivalent  eigenmode 
amplitude  will  be  unaffected  by  either  the  initial  curvature  or  the  longi¬ 
tudinal  load. 

The  spine  can  assume  various  equilibrium  positions.  Each  one 
will  produce  a  different  shape  of  ics  initial  curvature.  This  results  in 
different  sets  of  gQn  values.  The  amplitudes  of  the  various  eigenmodes 
and  locations  of  maximum  bending  stress  will  therefore  be  dependent  on 
the  attitude  of  the  subject. 

In  particular,  the  interpretation  of  test  results  with  animal  sub¬ 
jects  should  take  differences  in  curvature  between  man  and  animal  into 
consideration. 

C.  Dynamic  Instability 

When  a  time  variable  force  ia  added  to  the  longitudinal  load,  then 

P(t)*P+P0sinwt  .  (29) 

The  solution  that  was  just  discussed  does  not  apply.  The  differential 
equation  for  the  time  distribution  is  now 


f 

(P+P0  slnot)  -2-,  . 

Ml  Z  (30) 
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This  equation  represents  a  Mathieu  equation  except  for  the  addition  of  a 
forcing  function.  Due  to  the  multiplicative  connection  of  T  and  its  co¬ 
efficient,  the  solution  contains  also  harmonics  that  are  of  different  fre¬ 
quency  than  the  parameter  variation.  Which  one  of  these  harmonics 
dominates  the  solution  depends  on  the  relation  of  the  parameter  fre¬ 
quency  to  the  natural  frequency.  Unbounded  solutions  can  be  expected 

3 

at  and  around  frequencies  that,  according  to  Bolotin  ,  are 


,  k  =  1,  2,  3  . . . 


(31) 


The  solution  i3  of  frequency  «  if  k  is  even  and  of  frequency  «/2  if  k  is 
odd.  Therefore  there  are  frequencies  of  parametric  resonance  in  addi¬ 
tion  to  the  regular  resonant  frequencies.  These  frequencies  are  load- 
dependent.  The  relationship  between  the  first  resonant  frequency  and 
the  longitudinal  load,  equation  (20),  is  represented  by  the  solid  para¬ 
bola  (Figure  4)  that  intersects  the  abscissa  at  o>/q  =  1.  The  other 
parabolas  represent  locations  of  possible  parametric  resonance. 

Instability  is  also  possible  in  areas  around  these  parabolas.  Their 
width  depends  on  the  amplitude  of  the  load  variation.  These  regions 
of  instability  are  indicated  by  the  shaded  areas  in  Figure  4.  Most 
dangerous,  is  the  first  region  of  instability  that  is  represented  by  the 
first  region  on  the  right. 

Approximate  equations  for  the  boundaries  between  the  stable  and 

3 

unstable  regions  are  given  by  Bolotin  .  They  are: 

1st  Region  of  Instability: 


*>  *2  SX 


(32) 
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2nd  Region  of  Instability: 


*fl  Vu-IThhTT* 


pe” ”  3 


w.a^/n  -  4,(,“2'2) 

t 

3rd  Region  of  Instability: 


(33) 


(34) 


*  * aYp-p  j  *  PE  Euler  buckling  Load  . 
6 


Parametric  resonances  were  occasionally  observed  when  the  trans¬ 
verse  motion  of  the  spine  of  a  human  subject  under  vertical  vibrations 
was  observed.  The  difference  of  impedance  curves  as  obtained  through 
vibration  and  impact  tests  could  perhaps  be  explained  through  para¬ 
metric  excitation. 

Also  the  capability  of  the  spine  to  transfer  energy  at  various  fre¬ 
quencies  may  be  of  some  consequence,  in  particular  with  respect  to 
subsystems  that  are  coupled  to  the  spine.  Conversely,  coupled  sub¬ 
systems  may  distort  parametric  resonances. 

UNRESTRAINED  SPINE 

The  model  will  now  be  extended  to  include  the  thorax  without  ex¬ 
ternal  restraints.  The  boundary  conditions  at  the  upper  end  of  the  rod 
are  determined  by  the  inertia  of  the  thoracic  mass  to  transverse  accel¬ 
erations  and  the  rotational  inertia  of  the  thorax  (Figure  1).  A  longitu¬ 
dinal  force  is  acting  at  the  lower  end  of  the  lumbar  spine.  The 
configuration  is  somewhat  idealized  because  whole -body  accelerations 
and  rotations  are  assumed  to  be  negligible.  The  equation  of  motion  is 
the  same  as  in  the  previous  case  and  the  treatment  will  follow  the 
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same  sequence  of  steps  as  before.  The  homogeneous  equation  of  motion, 
equation  (2),  (straight  rod)  will  be  solved  first  and  the  effects  of  the 
longitudinal  loading  will  be  discussed.  The  second  step  is  the  solution 
of  the  nonhomogeneous  equation  of  motion,  equation  (1),  and  a  discus¬ 
sion  of  the  results. 

A.  Straight  Rod 

Equation  (2)  applies  in  this  case  as  well  as  the  ordinary  differen¬ 
tial  equations  (4)  and  (5)  that  were  derived  from  it,  and  the  general 
solution  for  the  spatial  distribution,  equation  (6),  and  time  distribution, 
equ  tion  (7).  The  constants  Dmn  will  have  to  be  determined  by  the 
boundary  conditions  which  were  discussed  in  the  previous  paragraph, 
and  are  quantitatively  defined  by  the  following  equations: 

V(o)  *  v"  (o)*0  ,351 


) 

X*| 


«  — m( 


) 

x  *! 


in-  mass  of  thorax 


(36) 


e  ( 


d3v 

dxdt2 


) 

x*l 


) 

x«t 


0  —dynamic  moment  of  inertia 
of  thorax 


(37) 


The  characteristic  equation  as  obtained  through  the  usual  manipulation 
is 
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Figure  5.  Curves  of  Possible  Instability  for  the  Unrestrained 
Spine  Model. 
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(38) 


In  this  equation  are 
,  m«n2 

8  =  ~sr 


(39) 


X 


El 


(40) 


The  characteristic  equation  contains  the  natural  frequencies  as  well  as 
the  longitudinal  load  P  which  is  contained  in  r1  and  (see  equation  7). 
The  natural  frequencies  can  therefore  be  determined  as  a  function  of  the 
longitudinal  load.  This  was  done  by  iteration.  The  results  are  presented 
in  Figure  5  by  the  solid  curves.  Arbitrary  values  were  chosen  for  the 
parameters.  The  only  purpose  of  these  preliminary  numerical  computa¬ 
tions  was  to  obtain  a  graphic  picture  of  the  relationship  between  load  and 
natural  frequency.  The  numerical  values  for  these  computations,  such  as 
the  modulus  of  elasticity,  static  moment  of  inertia,  etc.  were  rough  es¬ 
timates  of  the  various  spine  materials. 

The  eigenfunctions  for  this  set  of  boundary  conditions  are 

f*?eo»,aJI+8'i",a! 

- si.hr^  •  (41) 

fjjcothr^  —  fcsinbr^ 


These  eigenfunctions  depend  also  on  the  longitudinal  load.  The  first 
eigenfunction  at  three  different  longitudinal  loads  is  presented  in 
Figure  6.  The  numerical  values  for  the  various  parameters  were 
chosen  to  represent  a  human  subject.  The  length  f  is  30. 48  cm  which 
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2^ 

represents  the  length  of  the  average  lumbar  spine.  The  lowest  of  the 
longitudinal  loads,  Pj ,  is  27  lbs.  or  12.  27  kp,  which  is  about  the  weight 
of  the  average  human  thorax.  The  load  P2  is  three  times  that  of  P!  and 
the  load  P3  is  ten  times  that  of  P3.  It  is  very  obvious  that  these  eigen¬ 
functions  depend  very  much  on  the  longitudinal  load.  The  curve  for  the 
load  P2  is  in  shape  similar  to  the  transverse  displacement  distribution 
along  the  lumbar  spine  of  sitting  human  subjects  under  vertical  vibra¬ 
tion  (Figure  2). 

The  bending  moment  along  the  spine  in  one  of  its  eigenmodes  is  di¬ 
rectly  proportional  to  the  second  derivative  of  the  eigenfunction.  These  de¬ 
rivatives  are  plotted  in  Figure  7  for  the  same  three  longitudinal  loads.  It 
is  very  obvious  that  these  derivatives  depend  on  the  longitudinal  load.  The 
bending  load  is  rather  evenly  distributed  over  the  length  of  the  lumbar 
spine  at  small  longitudinal  loads.  Peaks  of  bending  loads  appear  as  the 
longitudinal  load  increases.  The  greatest  peak  occurs  at  the  upper  end 
and  a  smaller  one  at  the  lower  section  of  the  lumbar  spine.  The  peaks 
increase  with  the  longitudinal  load  and  the  lower  peak  moves  farther  down 
the  spine. 

Statistical  curves  have  been  published  that  show  the  incidence  of 
vertebral  injury  along  the  lumbar  and  thoracic  spine.  These  curves 

4 

differ  somewhat  from  one  author  to  the  other.  Hirsch  and  Nachemson 

arrive  at  a  distribution  along  the  lumbar  spine  that  is  fairly  even. 

5 

Moffatt  and  Howard  present  one  with  an  incidence  of  injury  that  is  low 

at  the  lower  end  of  the  lumbar  spine,  increases  slowly  from  there  on  up, 

and  increases  rapidly  around  the  11th  and  12th  thoracic  vertebra. 

6 

Another  curve  published  by  Higgins  et  al.  is  similar  to  that  of  Moffatt 
except  for  no  injuries  around  the  4th  lumbar  vertebra  and  some  injury 
of  the  coccyx. 

It  seems  that  the  differences  of  these  curves  may  be  due  to  differ¬ 
ences  in  the  magnitude  of  the  longitudinal  loads. 
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After  comparison  with  the  second  derivative  of  the  eigenfunction, 
one  would  expect  the  incidence  curve  of  Hirsch  and  Nachemson  to  be 
obtained  from  cases  of  relatively  small  longitudinal  loading.  The  sec¬ 
ond  derivative  of  the  eigenfunction  crosses  the  zero  line  for  two  longitu¬ 
dinal  loads.  This  seems  to  reflect  the  zero  incidence  in  Higgins'  curve 
around  the  4th  lumbar  vertebra. 


B.  Curved  Rod 

In  the  previous  discussion,  we  did  not  consider  the  right  hand  side 
of  equation  (1).  The  nonhomogeneous  equation  cannot  be  solved  exactly 
by  making  a  product  assumption  because  it  is  not  self  adjoint.  An  ap¬ 
proximate  solution  with  such  an  assumption  can  be  obtained.  In  this 
case,  the  Galerkin  method  will  be  applied. 


El 


The  equation  of  motion  is  rearranged  to  assume  the  form 

. 


d4v 


dt* 


+  P(t> 


dx‘ 


M  -4- P(t)  ■  L  (»)•  0  .  (42, 

at2  d*2 


A  product  assumption  of  the  following  nature  is  made 


v(x,t) 


M 


EV„  U)  Tft(t) 


(43) 


where  the  VR(x)  are  the  previously  determined  eigenfunctions.  This  is 
substituted  into  equation  (42)  which  results  in 

L(v)*«U.th  (44 

<(x,  t)  is  zero  if  such  an  assumption  can  satisfy  equation  (42)  exactly, 
and  it  is  different  from  zero  if  it  cannot  do  that.  In  this  case  the  error 
of  an  approximate  solution  with  this  assumption  is  minimized  by  im¬ 
posing  the  conditions 


M3 


* 

/l(v) VnU)  dx« 


n  =  1,  . . . .  N 


This  requires  the  error  «(x,  t)  along  the  length  of  the  rod  weighted  by 
the  eigenfunctions  and  integrated  over  x  to  be  zero.  In  other  words, 
equation  (45)  requires  orthogonality  between  the  residual  and  the  eigen¬ 
functions. 

The  variable  x  is  eliminated  through  integration  and  a  system  of 
coupled  ordinary  differential  equations  with  t  as  the  variable  evolves. 
These  equations  would  not  be  coupled  if  the  eigenfunctions  of  the  ad¬ 
joint  problem  would  have  been  used  for  weighting  functions.  The  coupled 
equations  are  in  matrix  notation 

{T„}  +-jr[ani]"'(EI  [bn,]  +  Pt,>[cni]  ){T„>. 


»nj  */ 

VflV) 

0 

I 

.* 

*  t 

d  *0 

0 

•77* 

The  solution  of  this  equation  consists  of  the  solution  for  the  homo¬ 
geneous  equation  plus  a  particular  solution  of  the  nonhomogencou* 
equation.  An  assumption  of  the  type 


{T0  *  (Tio) 


for  the  homogeneous  equations,  and  the  usual  manipulations,  leads  to 
a  solution  of  the  type 


(49) 


N 

T»  * 

k=  1 

The  are  the  amplitude  ratios. 

The  right  side  of  equation  (46)  is  constant.  Therefore,  a  constant 
can  be  determined  as  a  particular  solution.  This  solution  is 


{bJ  >  -  P<«)  [Hi-  [.yl-  {*,} 

(50) 

l«]  =  f  bijl  ''  (Elbijl  4  Pit)  [Cy]) 

(51) 

The  complete  solution  is  then 

»<«••)  ■  l\  l  «jk  „  *i»  (Uk,  4  *k)  *  Bn) .  (5M 

n=l  k=l 

This  can  be  rearranged  and  expanded  to  assume  the  following  form: 
V(x,t)  s  q,  +  Vjpa  «•  ..  ain(w|t  +  4*,) 

+  \  r^Hin  ♦  v#m  ♦  **  vnHnnJ  *in  («nl  *  V 

•* 

*  X  v  &  *  {53) 

n  n 

nal 

The  eigenfunctions  of  the  curved  rod  under  the  specified  boundary  con¬ 
ditions  can  be  expressed  as  series  of  the  eigenfunctions  of  the  straight 
rod  weighted  by*  the  amplitude  ratios  as  obtained  during  the  solution 
for  the  time  distribution. 

The  q  have  to  be  determined  using  the  initial  conditions.  Various 
n 

methods  can  be  applied  such  as  subdomain,  Galerkin.  or  a  collocation 
technique.  In  any  case,  the  q  will  finally  be  dependent  on  the  V  ,  that 
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is  the  eigenfunctions,  and  on  the  B  which  are  determined  by  the  initial 

n 

curvature.  Therefore  the  ^  .plitudes  of  the  free  motions  are  affected 
by  both  the  initial  curvature  and  the  longitudinal  load  which  is  also  part 
of  the  Bn's  and  the  eigenfunctions.  The  initial  curvature  of  the  spine 
changes  between  different  attitudes  of  a  subject.  The  difference  in  re¬ 
sponse  between  attitudes  can  be  accounted  for  through  the  initial 
curvature. 

C.  Dynamic  Instability 

The  solution  for  the  spatial  and  time  distribution,  equation  (49), 
applies  only  if  P(t)  is  a  constant.  The  differential  equation  for  the  time 
distribution,  equation  (46),  will  have  to  be  solved  by  procedures  such  as 
Hill's  method.  Of  particular  interest  are  areas  of  possible  dynamic  in¬ 
stability.  These  ha.  e  r.ot  been  determined  for  this  particular  case. 
However,  these  areas  are  usually  around  the  lines  of  the  natural  fre¬ 
quencies  in  Figure  5  and  also  around  curves  of  twice  the  natural  fre¬ 
quency,  fractions  of  it,  and  at  frequencies  that  are  the  sums  and 
differences  of  the  natural  frequencies.  Some  of  these  conditions  of 
parametric  instability  are  indicated  by  the  broken  lines  in  Figure  5. 

CONCLUSIONS 

It  was  the  purpose  of  this  discussion  to  investigate  the  possibilities 
of  accounting  quantitatively  for  the  effect  of  longitudinal  loads  on  the 
natural  frequencies  of  transverse  motions  of  the  lumbar  spine,  as  well 
as  for  the  effects  of  longitudinal  loading  and  initial  curvature  on  the 
magnitude  and  distribution  of  bending  along  the  lumbar  spine. 

Longitudinal  loads,  initial  curvature  and  boundary  conditions  im¬ 
posed  on  the  lumbar  spine  seem  to  have  significant  effects  on  the  mag¬ 
nitude  and  distribution  of  bending  in  a  dynamic  environment.  The  type 
of  bending  distribution  curves  that  have  been  obtained  seem  to  agree 
with  statistical  curves  of  injury  incidence  along  the  lumbar  spine.  The 
transverse  displacement  distribution  along  the  lumbar  spine  under 
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longitudinal  loads  agrees  with  experimental  data.  A  more  detailed 
numerical  adjustment  of  the  model  to  experimental  data  will  be  attempted. 
In  particular,  major  subsystems  that  are  coupled  to  the  spine  will  have 
to  be  considered  for  this  purpose.  Because  subsystems  have  an  effect 
on  the  boundary  conditions  and  thus  on  the  response  of  the  spine,  it 
seems  they  have  to  be  included  because  no  experimental  data  are  avail¬ 
able  that  exclude  coupling  effects. 

The  natural  frequency  of  transverse  spine  motions  under  longitudinal 
loads  is  load-dependent.  The  natural  frequency  decreases  in  general 
with  increasing  compressive  loading  and  increases  with  increasing  ten¬ 
sile  loading.  The  type  of  natural  frequency  -  longitudinal  load  relation¬ 
ship  depends  to  a  great  degree  on  the  boundary  conditions. 

Parametric  transverse  resonances  and  dynamic  instabilities  can  be 
expected  under  a  periodic  longitudinal  load. 
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